Weak ferromagnetism in Mn nanochains on the CuN surface 
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We investigate electronic and magnetic structures of the Mn chains supported on the CuN sur- 
face using first-principle LSDA and LDA-I-U calculations. The isotropic exchange integrals and 
anisotropic Dzyaloshinskii-Moriya interactions between Mn atoms are calculated using Green func- 
tion formalism. It is shown that the account of lattice relaxation and on-site Coulomb interaction 
are important for accurate description of magnetic properties of the investigated nanosystems. We 
predict a weak ferromagnetism phenomenon in the Mn antiferromagnetic nanochains on the CuN 
surface. The value of a net magnetic moment and direction of spin canting are calculated. We show 
that some experimental features may be explained using anisotropic exchange interactions. 

PACS numbers: 73.22.-f, 75.30.Et 
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I. INTRODUCTION 

The investigation of magnetic nanomaterials is impor- 
tant part of nanoscience and exerts influence on progress 
in different sectors of technology (such as medicine de- 
vices for therapy and diagnostics,^'^ magnetic data stor- 
age systems;^ etc). In the presence of intrinsic mag- 
netic moment scientists can change a physical proper- 
ties of nanomaterials by applying external magnetic field. 
The technological applications of magnetic nanomateri- 
als should base on an accurate control of the coupling 
between individual spins. 

Recently, Hirjibehedin et al.'^ have reported the fabri- 
cation of Mn nanoparticles in the form of linear chains 
that display truly collective quantum behaviour. Using 
local spin-excitation spectroscopy technique, based on in- 
elastic scanning tunneling microscopy (STM), they were 
able to show how the quantum properties of this sys- 
tem depend on the number of atoms involved. They 
demonstrated an innovative method to measure and con- 
trol these magnetic interactions. The experimental spec- 
trum was analyzed using the simplest form of Heisenberg 
model with exchange interaction only between nearest 
neighbors. However, there are a number of experimental 
results which cannot be explained by authors of Ref|4| 
using the Heisenberg model: (i) zero-field splitting which 
grows in energy with increasing chain length, (ii) the dif- 
ferent zero- field energy of the m = ±1 and m — excited 
states and (iii) asymmetry of the spectra with respect to 
voltage polarity. 

Jones and Lin^ have applied GGA-I-U approach to de- 
scribe the electronic and magnetic structures of single 
and pair of Mn atoms on the CuN(lOO) surface. The 
performed spin-density analysis shows that Mn atoms on 
such surface preserve their atomic spins S — ^. This 
result agrees with STM measurement."*^ Electron-density 
change and surface relaxation due to Mn atoms are also 
analyzed in Ref!^. 

The combination of experimental STM approach and 



theoretical ab-initio methods^ has been used in order to 
describe the large magnetic anisotropy of individual Mn 
and Fe atoms on the CuN surface. The authors of the 
paper^ have provided the detailed phcnomcnological pic- 
ture of magnetic anisotropy and concluded that in case of 
manganese system the easy axis is oriented out-of-plane. 

In this paper we show that the local distortion of the 
system results in a superexchange interaction between 
Mn atoms through N atoms. Isotropic exchange interac- 
tions are calculated using Green functions approach and 
total energies difference method. Using full diagonaliza- 
tion of Heisenberg Hamiltonian with calculated isotropic 
exchange integrals we estimate the energies of first mag- 
netic excitations. The results are in good agreement with 
experimental data. 

In the previous theoretical investigations a non- 
coUinear magnetic ground state for nanostructures on 
non-magnetic surfaces^i^ and magnetio^ were proposed. 
For instance, the results for the non-collinear triangu- 
lar compact trimer of Cr on the Au(lll)^ predict that 
the angle between each pair of moments equals to 120° 
and the total spin moment is zero. Therefore, the non- 
collinearity is result of frustration of magnetic interac- 
tions. 

In the paper— authors have investigated different ge- 
ometries of Fe, Mn and Cr atoms on the Cu(lll) sur- 
face. The Fe clusters were found to be ferromagnetically 
ordered. Whereas for the Mn and Cr clusters an antifer- 
romagnetic exchange interactions between nearest neigh- 
bours have been found. The antiferromagnetic couplings 
produce either coUinear or non-collinear magnetic struc- 
tures due to frustration of cluster geometry. 

An interesting results for the trimer and tetramer con- 
figurations of Mn and Cr atoms on the magnetic Ni(lll) 
surface were obtained in the paper One should stressed 
that there are two types of a magnetic frustration: (i) 
frustration within adcluster and (ii) frustration arising 
from competing magnetic interactions between the ad- 
clusters and the surface atoms. 
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Thus, one can conclude that the only known source 
for spin non-coUinearity of magnetic clusters on non- 
magnetic 3d surface is geometrical frustration which re- 
sults in magnetic frustration. In this paper we pro- 
pose a new source of spin non-collinearity for nanosys- 
tems on a surface. According to our calculations a lo- 
cal distortion between Mn atoms in the nanochain re- 
sults in the strong Dzyaloshinskii-Moriya (DM) interac- 
tion. An important role plays the displacement of N atom 
from the surface. Based on first-principle calculations of 
the Dzyaloshinskii-Moriya interactions between magnetic 
moments we point out that the Mn nanochains on the 
CuN demonstrate a weak ferromagnetism phenomenon. 
We have estimated the value of a net magnetic moment 
and direction of spin canting. These results are also con- 
firmed by direct LDA-I-U+SO calculations. 

The paper is organized as follows. In Section II we 
describe the methods of the investigation. In Section III 
A and III B we present the results of LSDA and LDA+U 
calculations, respectively. The analysis and comparison 
of obtained exchange interactions with experimental data 
are presented in Section III C. Section IV is devoted 
to the analysis of zero-field energy splitting observed in 
STM experiment and in section V we briefly summarize 
our results. 



II. METHODS OF INVESTIGATION 

A. DFT calculation details 

We have used two complementary approaches for in- 
vestigations of an electronic and magnetic properties of 
Mn nanochains on the CuN surface. 

(i) First-principles total-energy and force calculations 
were carried out using the projected augmented-wave 
(PAW) methodic as implemented in the Vienna ab initio 
simulation package (VASP)?- ^-" ' ' — Exchange and correla- 
tion effects have been taken into account using LSDA and 
LDA-I-Ui^ approaches. In all cases under investigation 
we used an energy cutoff of 400 eV in the plane-wave 
basis construction and the energy convergence criteria 
of 10~* eV. The atomic positions of considered systems 
were relaxed with residual forces less than 0.01 eV/A. For 
the Brillouin zone integration, a (4x4x1) Monkhorst-Pack 
mesb'^'* and Gaussian-smearing approach with cr = 0.2 eV 
were used. 

To simulate structure of the unit cell we have used a 
supercell approach. Structure of the supercell has con- 
sisted of two-layer (2 x (n-t-1)) Cu(lOO) surface, N atoms 
embedded into upper Cu-layer, Mn„-chains placed on the 
top of the CuN-surface and vacuum region of 10 A. Lat- 
tice constant for Cu was chosen to be 3.63 A, which gives 
a minimal value of the total energy in calculation of the 
bulk fee Cu. Lower layer of Cu has been fixed under 
relaxation. 

(ii) We have also used the Tight Binding Linear-Mufin- 
Tin-Orbital Atomic Sphere Approximation (LMTO) 



methodic in terms of the conventional local den- 
sity approximation taking into account the on-site 
Coulomb interaction LDA-I-U and spin-orbit coupling 
LDA-|-U-|-SOji^iii In this type of calculations we have 
used the relaxed structures obtained by PAW approach. 
The radii of atomic spheres were r(Mn) = 1.137 A, 
r(Cu) = 1.322 A and r(N)=0.793 A. In order to fiU the 
empty space of the unit cell required number of empty 
spheres were added. 

B. Spin Hamiltonian approach 

The main aim of our investigation is first-principle de- 
termination of parameters of the following spin Hamilto- 
nian: 

H = Hueis + Hdm, (1) 
where HeisenbergiS energy term is 

Hueis — ^^Jij^i^ji (2) 
i<j 

and Dzyaloshinskii-Moriyai^ energy term is 

iJiJM = ^Aj [^, X 4]. (3) 

In order to calculate the isotropic exchange interac- 
tions 3ij, in Eq.(2) between magnetic moments of Mn 
atoms we have used two different approaches, (i) The to- 
tal energy difference method on the basis of PAW results. 
The main idea of this method is that the isotropic ex- 
change interaction defines through the energy differences 
between different magnetic configurations. For instance, 
in case of Mn-dimer, spin Hamiltonian of the system can 
be written in the following form: 

H^JSi-S2. (4) 

The corresponding total energies of the ferromagnetic 
and antiferromagnetic configurations of two classical 
spins are given by 

Efm = JS^ (5) 

and 

Eafm = -JS^ (6) 

Therefore, the exchange interaction J is expressed in the 
following form: 

J _ EpM — Eafm ,„x 

(ii) From the other hand, based on LMTO results 
one can calculate the isotropic exchange integrals and 
Dzyaloshinskii-Moriya interactions between magnetic 
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moments of Mn atoms (S=|) using the local force theo- 
rem and Green functions formalism^SiiSii^^ 
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where m (m , m , m ) is magnetic quantum number 
and the on-site potential A™™' = 
Green function is calculated in the following way: 
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FIG. 1: The projection on yz-plane of relaxed structure 
of Mn-dimer supported on the CuN obtained using LSDA 
approach. Yellow, blue and red spheres correspond to Cu, 
N and Mn atoms, respectively. Green arrows correspond to 
direction of atomic maenetic moments. 



Here c™" is a component of the n-th eigenstate, E" is 
the corresponding eigenvalue and k is quasimomentum 
in the first Brillouin Zone. 

In turn the Dzyaloshinskii-Moriya interaction, Eq.(3) 
can be calculated through the account of spin-orbit cou- 
pling in the second variation of total energy of the system 
over the small deviations of magnetic moments from the 
coUinear ground statei^i^ 
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where — XkLS and Xk is spin-orbit coupling con- 
stant for site k. Here we present only z component of 
Dzyaloshinskii-Moriya vector, x and y components can 
be obtained from the z ones by rotation of the coordinate 
system. 
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FIG. 2: Total and partial density of states obtained using 
LSDA calculations, (a) Total DOS, (b) Projected DOS onto 
Cu atoms, (c) Projected DOS onto N atoms, (d) Projected 
DOS onto Mn atom. Fermi level corresponds to eV. 



III. RESULTS 
A. LSDA results 

We have performed LSDA calculation of antifcrromag- 
netic Mn-dimer supported on the CuN surface. Fig[T] 
shows the relaxed structure of the Mn-dimer obtained 
within LSDA method using PAW approach. The infor- 
mation about structure of relaxed system is presented in 
Table I. One can see that Mn-N-Mn bond angle of 171° 
is close to 180° and corresponds to the maximum of su- 
perexchange interaction between 3d atoms. 

The calculated total and partial densities of states ob- 
tained using PAW approach are presented in Fig. [21 The 
valence band contains low and high energy parts which 
are separated by 6 eV. The low energy states located 
around -15 eV are mainly N ones. In the region from -7 
eV to 5 eV all states are highly mixed. 



The calculated values of magnetic moments of Mn 
atoms within LMTO and PAW approaches are 3.70 y^B 
and 3.35 /is, respectively. These values are smaller than 
experimentally observed spin |. Moreover, the calcu- 
lated exchange parameter within Green functions ap- 
proach, Eq.(8) is 20.4 meV, whereas total energies dif- 
ference method value, Eq.(7) is 24.8 meV. These values 
at least three times larger than that experimentally ob- 
served. Thus, one can see that the main electronic and 
magnetic properties of the investigated nanosystem can- 
not be correctly reproduced within the LSDA approach. 



B. LDA+U results 

The results of the previous section have demonstrated 
drawbacks of the LSDA approach to describe the mag- 




FIG. 3: The yz projections of relaxed structures of Mn chains 
of CuN surface for n=l ^ 4. Yellow, blue and red spheres 
correspond to Cu, N and Mn atoms, respectively. Green ar- 
rows correspond to direction of atomic magnetic moments of 
manganese. 



netic properties of the Mn dimer on the CuN surface. 
It is well known problem of local density approximation 
in respect to transition metal compounds. To overcome 
this problem we have used the LDA+U approach with 
on-site Coulomb and on-site Hunds interaction parame- 
ters of U=6.0 eV and J^r=0.9 eV, respectively. These 
values are in good agreement with recent first-principle 
estimations performed in the worki^ 

The relaxed structures of the Mn-chains of different 
lengths (n=l 4) obtained within PAW calculations are 
shown in Fig [31 The structural information is presented 
in Table HI The obtained structures have some inter- 
esting geometrical features. Let us analyze the differ- 
ence between the clean CuN-substrate and the substrate 
with Mn adatoms on the top. Presence of Mn atoms 
causes some rearrangement of upper layer atoms of CuN- 
substrate. Generally, this rearrangement concerns the N 
atoms. In contrast to clean-CuN surface, the N atoms of 
the system with Mn nanochains are significantly shifted 
from the first layer plane. This fact agrees with results of 
recent GGA calculations.^ The calculated angle of Mn-N- 
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FIG. 4: Total and partial density of states of Mn-dimer on 
CuN surface obtained using LDA-(-U calculations, (a) Total 
DOS, (b) Projected DOS onto Cu atoms, (c) Projected DOS 
onto N atoms, (d) Projected DOS onto Mn atom. Fermi level 
corresponds to eV. 

Mn bond within LDA+U approximation equals to 143°. 
The N atoms at the edges of chains also have some dis- 
placement from the plane in z-direction, but to a smaller 
extent than N atoms situated inside the chain. 

From a geometrical point of view the important differ- 
ence between results of LSDA and LDA-I-U approaches is 
the angle of Mn-N-Mn bond. Let us analyze this fact on 
the level of hopping integral. For simplicity, we assume 
that there is the only strong hopping between orbitals of 
two 3d atoms. The hopping integral tij is proportional 
to cos a, where a is angle of metal-ligand-metal bond. 
Therefore, within LDA+U approach the hopping integral 
is strongly suppressed due to local Coulomb correlations. 
In turn the isotropic exchange interaction between Mn 
atoms in the atomic limit of Hubbard model can be ex- 
pressed as Jy = here U is on-site Coulomb integral. 
It is clear that on-site Coulomb interaction is another 
source of suppression of the isotropic exchange interac- 
tion in LDA+U in comparison with LSDA approach. 

The total and partial densities of states of the dimer 
and trimer systems obtained using LDA+U approxima- 
tion are presented in Fig|4] and [5l respectively. The cal- 
culated value of magnetic moment are listed in Table 
im One can see that the values of magnetic moments of 
Mn atoms within the chains varies from 4.34 (middle 
atoms) to 4.47 (edge atoms) and now is much closer 
to experimental values of 5 = |. In the case of LMTO 
results this difference is much smaller. 



C. Isotropic exchange interaction 

The next step of our investigation is determination of 
Heisenberg exchange interaction parameters in Eq.(2). 
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FIG. 5: Total and partial density of states of Mn-trimer chain 
on CuN surface obtained using LDA+U calculations, (a) To- 
tal DOS, (b) Projected DOS onto Cu atoms, (c) Projected 
DOS onto N atoms, (d) Projected DOS onto central Mn atom, 
(e) Projected DOS onto edge Mn atom. 

TABLE I: Structural information about alignment Mn and N 
atoms on the Cu(lOO) surface as result of relaxation in cases 
of LSDA and LDA+U approaches. The values are distance 
(in A) between Mn atoms, ^-coordinates of Mn atoms and N 
atoms at the center and edge of chain, respectively, a is angle 
of Mn-N-Mn bond. Zero level of z coordinate corresponds to 
the lower layer of Cu atoms. All values are given in A. 
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3.58 


3.66 


3.49 


2.66 


171° 
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3.78 
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TABLE IL Values of magnetic moments of Mn atoms (in /is) 
calculated using LMTO (PAW) method. 

n Medge M 

center 

~ '- 4.46 (4.48) 

2 4.45 (4.44) 

3 4.47 (4.47) 4.45 (4.34) 

4 4.47 (4.47) 4.40 (4.35) 



The magnetic couplings between Mn atoms calculated 
using Green functions, Eq.(8) and the total energy differ- 
ence method, Eq.(7) are presented in Table Hill One can 
see that the calculated values of the dimer interaction are 
in good agreement with experimental value of 6.4 meV. 
The value of isotropic exchange integral between nearest 
Mn atoms in trimer is smaller than that in dimer system. 
There is also small ferromagnetic coupling between edge 
Mn atoms. 

In the case of quatromer system there is the difference 
between exchange integrals J12 and J23. The coupling 
at the center of the chain is smaller than the coupling 
at the edge. There is interesting dependence of nearest- 
neighbour exchange interaction according to the chain 



length. This tendency probably corresponds to oscilla- 
tions of exchange interaction parameter depending on 
chain length. It is important to investigate the mecha- 
nism of such strong oscillations (see Mns results in Table 
imp of nearest neighbour exchange interactions in Mn„- 
nanochains. Such analysis is left for future investigation. 

Using full diagonalization procedure of ALPS 
library^i^ we have calculated the spin excitation spec- 
tra of investigated systems. The energies of first excited 
states are presented in Table HVl Despite of the fact that 
the J12 of trimer has smaller value than in dimer case, 
weak ferromagnetic J13 interaction between edge atoms 
compensates this difference and gives us opportunity 
to reproduce experimentally observed excitation energy 
with reasonable accuracy. 

TABLE III: Values of exchange interactions Jij (in meV) be- 
tween magnetic moments of Mn atoms calculated using Green 
function method (TB-LMTO-ASA) method. Values obtained 
using total energies difference method (PAW) approach are 
given in parenthesis. 

n J12 J13 J23 J34 J24 

2 7.0 (6.0) _ _ _ _ 

3 4.0 (4.2) -0.09 (-0.09) 4.0 (4.2) 

4 5.6 (5.2) -0.07 (-0.04) 2.4 (4.1) 5.6 (5.2) -0.07 (-0.04) 



For quatromer system our results are in excellent 
agreement with experimental spectrum. 

TABLE IV: Energies of first excited states (in meV) of Heisen- 
berg model obtained using ALPS code for the chain systems 
of different length. 



n 




^LMTO 


TpcaLc 
t^PAW 


2 


6.4 


7.0 


6.0 


3 


16.0 


10.5 


10.9 


4 


2.9 


3.0 


2.6 



D. Dzyaloshinskii-Moriya interaction 

From the crystal symmetry point of view there is no 
inversion center at the point bisecting the straight line be- 
tween Mn atoms of investigated nanosystems. Therefore, 
in according with Moriya's rules^^ a DM coupling exists. 
First, let us perform a simple geometrical analysis of the 
symmetry of Dzyaloshinskii-Moriya vector. There are 
two sources of inversion symmetry breaking in the inves- 
tigated systems, (i) The first one is the substrate surface. 
Based on the fact that previous investigations of metal- 
lic nanochains on nonmagnetic 3d surface has no sign 
of non-coUinearity, one can conclude that surface gives 
negligible small contribution to anisotropic exchange in- 
teraction, (ii) More importantly, the second source of 
inversion symmetry breaking is vertical displacement of 
N atom. The Dzyaloshinskii-Moriya vector, _Di2 between 



Mn atoms is proportional to [f x R12] where R12 is a 
unit vector along the line connecting the magnetic ions 
and r is the shift of the ligand atom from this line (Fig. 
[S]). One can see that in our case r and R12 have z and 
y components, respectively. Therefore, the direction of 
Dzyaloshinskii-Moriya vector is x axis. 



R12 



Mn 



r 



Mn 



FIG. 6: Schematic representation of displacement of the lig- 
and atom. 

The calculated anisotropic exchange interactions, 
Eq.(lO) are presented in Table IVl For all systems un- 
der consideration the Dzyaloshinskii-Moriya vector lies 
along X axis (perpendicular to the Mn chain and paral- 
lel to the CuN surface). Therefore, if all spins lie in the 
yz plane, the canting exists and there is weak ferromag- 
netism in the system. The ratio between Dzyaloshinskii- 
Moriya and isotropic exchange interactions, -^^^^=0.002 
is the same order of magnitude as in case of well known 
antifcrromagnets Fe203 and La2Cu04 with weak ferro- 
magnetism. 



TABLE V: Values of x component of Dzyaloshinskii-Moriya 
interactions Dij (in meV) between magnetic moments of Mn 
atoms calculated using Green functions method (Eq.(lO)) 
within LDA-I-U approach. 



n 


Df2 


Di3 D23 


2 


0.014 




3 


0.018 


0.000 0.018 


4 


0.024 


-0.006 0.030 



We have minimized the classical spin Hamiltonian 
(Eq.(l)) with first-principles exchange parameters in re- 
spect to the angle between different spins in the chain. 
On this basis we have defined the values of canting angles 
and net magnetic moments of the Mn nanochains. These 
results are presented in Table |Vl] and Fig. [T] One can 
see that the net magnetic moment increases with length 
of the nanochain. 

In order to test the reliability of the weak ferromag- 
netism results we have performed the LDA+U+SO cal- 
culations. For dimer system the magnetic ground state 
is non-collinear and spins are along z axis with canting 
of 1.0° for LMTO and 1.6° for PAW. These results are in 
good agreement with previous GGA calculations^ where 
the easy axis of a system with single Mn atom on the 
CuN surface is z axis. The obtained values of canting an- 
gles are about one order larger than those obtained in our 




FIG. 7: Schematic representation of spin canting of dimer, 
trimer and quatromer systems. 



LDA-I-U calculations f Table IVH) . This observation can be 
addressed to underestimation of Dzyaloshinskii-Moriya 
interaction calculated by Green functions method.— 
Obtained canting angles correspond to the following 
Dzyaloshinskii-Moriya interactions Dfj — Jijtg{'K — 2a) 
= 0.24 meV (in LMTO) and D^^ = 0.34 meV (in 
PAW) . One can consider these results as manifestation of 
strong Dzyaloshinskii-Moriya interaction in investigated 
nanosystems. 



TABLE VI: Values of canting angles and full weak ferromag- 
netic moments of different Mn nanochains (in ^b) obtained 
using minimization of model spin Hamiltonian. 



n 


angle 


^'imodel 


2 


a=0.057° 


0.009 


3 


a=0.164° 


0.038 


4 


a=0.173° /3=0.466° 


0.045 



Clearly, the ultimate test of our results will to com- 
pare them with experiment. In the next section we will 
show that some experimentally observed features- can be 
explained using the anisotropic exchange interaction. 



IV. ZERO-FIELD SPLITTING 



In order to explain experimentally observed different 
zero-field energy of to = ±1 and m = excited states^ we 
use the quantum spin Hamiltonian with Dzyaloshinskii- 
Moriya interaction: 



H ^ JSiS2 + Di2[Si X ^2]. 



(11) 
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For simplicity, let us consider the case of S=^ and 

Di2 — {D'i2\ 0; 0). One can rewrite Eq.(ll) in the follow- 
ing form: 

H = J(^f 5f + SlSl + + DUSIS^^ - 5J5|)(12) 

The basis functions for this Hamiltonian can be written 
as follows, 

IT,T> U,i> IT,i> U,T>- (13) 

Using well known rules for spin operators 



splitting. We plan to investigate the effect of single-ion 
magnetic anisotropy energy on this splitting. 

TABLE VIII: The energy difference, A£; of m = and m = 
±1 excited states for different sets of calculated isotropic and 
anisotropic exchange interactions (in meV). 



^^IT>=2U> 



^1T>=^IT> 



5"U>=^IT> ^^U>=-^IT> ^1i>=-^U> 

one can define the matrix elements of this Hamiltonian 
presented in Table IVIII The eigenvalues of this matrix 



TABLE VII: Matrix elements of the Heisenberg Hamiltonian 
Eq.(12) 



T,T> U,i> I T,i> I J.,T> 



1 T,T> 


J 

4 





»DJ2 
4 


4 


1 i,i> 





J 

4 


iDf^ 
4 


4 


1 T,i> 






J 


J 


4 


4 


4 


2 


1 i,T> 


4 


4 


J 

2 


J 
4 



are the following: 





4D2 


2 







set 




AE 


LMTO (LDA-hU): J = 


7.0, 0.02 


< 10-* 


LMTO (LDA+U+SO): J 


= 7.0, D'' = 0.24 


0.02 


PAW (LDA-I-U-I-SO): J = 


= 6.0, = 0.34 


0.04 



V. CONCLUSION 

We have performed first-principle investigations of 
electronic and magnetic structures of the Mn nanochains 
supported on the CuN surface. Relaxation effects have 
taken into consideration. The calculated isotropic ex- 
change integrals are in good agreement with experimental 
data. We have also calculated the anisotropic exchange 
interactions in the system and predicted the antiferro- 
magnetic ground state with weak ferromagnetism. We 
stress that the main source of this phenomenon is local 
distortion which breaks the inversion symmetry between 
Mn atoms. It follows that the relaxation effects are im- 
portant for the system under consideration. The calcu- 
lated values of canting angles are larger than those for 
classical antiferromagnetics with weak ferromagnetism, 
Fe203 and La2Cu04. Using calculated anisotropic ex- 
change interactions we have explained the experimen- 
tally observed different zero-field energies of m — and 
m — ±1 states. 

Based on our results one can expect the weak ferromag- 
netism phenomenon in the similar surface nanosystems. 
For instance, we found this spin-orbit coupling effect in 
Co nanochain on the Ft surface. Such work is in progress. 



One can see that the energies of m = and m = ±1 
triplet states are different. Therefore, one can expect that 
anisotropic exchange interaction helps us explain similar 
difference in the experimental spectra for the Mn dimer 
on the CuN surface. 

Since in the case of S=5/2 the situation is more 
complicated, we have numerically calculated the exci- 
tation spectra using previously obtained isotropic and 
anisotropic exchange interactions by means of ALPS 
Iibraryi2ii2^ The final results are presented in Table IVIIII 
According to experiment the zero-field energies of ex- 
cited states are 5.96 ± 0.05 meV and 5.83 ± 0.05 meV 
for m = and m = ±1, respectively, and correspond to 
AE — 0.13 ± 0.05 meV. One can see that the theoreti- 
cally estimated values of AE (Table IVIIip give the cor- 
rect order of magnitude for experimental zero-field energy 
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